Layering instability in a confined suspension fiow 
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We have shown [J. Fluid Mech. 592, 447 (2007)] that swapping (reversing) trajectories in confined 
suspension flows prevent collisions between particles approaching each other in adjacent stream- 
lines. Here we demonstrate that by inducing layering this hydrodynamic mechanism changes the 
microstructure of suspensions in a conflned Couette flow. Layers occur either in the near-wall 
regions or span the whole channel width, depending on the strength of the swapping-trajectory 
effect. While our theory focuses on dilute suspensions, we postulate that this new hydrodynamic 
mechanism controls formation of a layered microstructure in such flows in a wide range of densities. 

PACS numbers: xxxx 



Confined particulate flows are important due to 
their applications in microfabrication, microfluidics and 
biotechnology. Recent studies have shown complex non- 
linear microstructural evolution [ll-13] , strikingly different 
than the behavior with no confining walls. 

Confinement affects the system behavior both by im- 
posing geometrical constraints on particle motion and by 
giving rise to purely hydrodynamic mechanisms. While 
geometry-driven phenomena have been extensively stud- 
ied (e.g., ordered structures in dense suspensions and 
emulsions tightly confined by bounding walls l^), hy- 
drodynamic mechanisms that are wall-induced have only 
recently been at the center of attention. 

So far the most thoroughly investigated hydrodynamic 
confinement phenomenon is the fluid backflow produced 
by particle motion. It has been studied in linear conduits 
[3, 0| , and in parallel- wall channels [1, [l3| • The backflow 
gives rise to the anomalous sign of the mutual diffusion 
coefficient for Brownian spheres 0, [i3K wave propagation 
in flow-driven drop or particle trains 'Sj'B,'?!], and unusual 
stability of a square particle lattice _5, 7]. 

Here we demonstrate that the evolution of suspen- 
sion microstructure also depends on the swapping (re- 
versing) trajectory effect Because this proposed 
conflnement-induced hydrodynamic mechanism prevents 
collisions of particles in adjacent streamlines in a con- 
fined Couette flow, a uniform suspension microstructure 
is destabilized, leading to the formation of particle layers 
parallel to the walls (cf.. Figs. [T] and [5]). 

Layered suspension microstructure that spontaneously 
arises in conflned flows, even at moderate particle con- 
centrations, has recently been observed experimentally 
[H and in computer simulations (also in movies Al~ 



A4 in IJ]). A key hydrodynamic mechanism responsible 



for this behavior is revealed by our present work. 

We consider a suspension of non-Brownian spherical 
particles in planar Couette flow of shear rate 7 [cf.. Fig. 
[ija)], under creeping- flow conditions. The particles have 
finite roughness, which is modeled by a very steep, short- 
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FIG. 1. Wall-induced layering, (a) Suspension geometry, (b) 
Swapping-trajectory (ST) effect: the flow scattered from the 
walls produces lift causing the reversal of particle motion, (c) 
Relative particle trajectories in two systems: unconflned (left) 
and conflned with an ST region (right), (d) Time evolution 
of the transverse position z in suspension: unconflned (left) 
and conflned (right). Time is normalized by the characteristic 
time between collisions r — l/{nojd'^). 



range repulsive potential of the range dr — d + e (where 
d is the hydrodynamic diameter) to mimic contacts be- 
tween roughness asperities. These direct particle con- 
tacts remove the Stokes-flow symmetry of binary colli- 
sions, producing finite transverse particle displacements 
[l5| . as illustrated in Fig.[T][c). 

Our analysis is focused on the dilute-suspension limit. 
The simulations in Figs. [T] and [2] are performed using the 
Boltzmann-Monte Carlo (BMC) method, in which the 
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FIG. 2. Suspension microstructure for particles of different 
rougiiness in a ciiannel of widtii H/d = 15 at time tjr = 2400. 
Density profiles n/no are shown vs transverse position zjd. 
Roughness range (a) e/d = 0.25; (b) e/d = 0.064. 



system dynamics is modeled as a sequence of uncorre- 

Since particle correlations are 



lated binary collisions [11 



important only at higher concentrations, this approach is 
appropriate in the low-concentration regime. We analyze 
suspensions at low concentrations no to emphasize phys- 
ical mechanisms involved in the formation of the layered 
microstructure. However, our direct numerical simula- 
tions (cf., movies B1-B4 in [l3|) demonstrate that these 
mechanisms are present also at higher concentrations. 

In the BMC method, the particle distribution is 
evolved by performing a sequence of uncorrelated col- 
lisional steps. In each step, we choose a random pair of 
particles, simulate their binary collision, and update the 
cross-streamline particle positions according to the post- 
coUisional displacements IJ]. When binary trajectories 



are accurately evaluated taking into account the hydro- 
dynamic interactions (HI) in the wall presence [l3, [l3| , 
we refer to this description as BMC-HI. To highlight the 
role of topological features of pair trajectories without 
incurring large numerical cost of evaluation of HI, we use 
two simplified collision models M2 and M3 (described be- 
low). Figures [U and [2] were obtained using BMC-HI, and 
Figs. [SHU correspond to models M2 and M3. 

Figure [Hd) demonstrates that a suspension confined 
between two parallel walls develops a well-defined layered 
structure (after about 20 collisions per particle), while 
the unconfined system remains uniform. Fig. [2] illustrates 
a strong dependence of the layered structure on the range 
e of non-hydrodynamic interparticle interactions. 

We argue that the observed layering behavior stems 
from the swapping-trajectory (ST) effect that causes ap- 
proaching particles to reverse their motion, and avoid col- 
lision. Such trajectories, which occur in the wall presence 
[cf., Fig. [Ijc)], prevent large collisional displacements for 
particles with a sufficiently small transverse offset, lead- 
ing, therefore, to stabilization of particle layers. 

The ST domain results from wall-mediated HI between 
the approaching particles JlIJ . The wall reflection of the 
perturbation flow produced by one of the particles [cf.. 
Fig. [Hb)] pushes the other particle across streamlines of 
the applied flow toward the fluid moving in the oppo- 



site direction, causing the reversal of the relative particle 
motion. Through this purely hydrodynamic mechanism 
confining walls influence the suspension microstructure. 

Population-balance equation - To elucidate the effect 
of swapping trajectories on the suspension dynamics we 
use the population-balance method, where the particle 
density n{z) (uniform in the flow and vorticity directions 
X and y) is described by a master equation that accounts 
for the effect of binary particle collisions on the suspen- 
sion motion. In the simplest two-dimensional (2D) ver- 
sion (with no y direction), the master equation reads 



dn {z, t) 



7/ [n(z + iA-iA')n(z 
n{z)n{z - A)] AdA, 



(1) 



where A and A' are the pre-collision and post-collision 
particle offsets. For simplicity, we assume that the inter- 
acting particles undergo symmetric transverse displace- 
ments. The flrst term of the integrand on the right-hand 
side of Eq. (P) corresponds to collisions pushing a parti- 
cle into the position z (moving it from z -\- ^A — ^A'), 
and the second term to collisions displacing a particle 
from the position z. If the initial and final offsets are the 
same, A = A', the flrst and second terms of the inte- 
grand cancel. In the full three-dimensional (3D) version, 
an additional integration with respect to the lateral offset 



Ay would be present [14 1 



Collision models M2 and M3 - Based on geometry of 
binary collisions we introduce a 2D collision model M2, 
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FIG. 3. (a) Pair collisions: relative trajectories for roughness 
range e/d = 9 • 10~^ (dashed lines) and e/d = 10""^ (gray 
solid), (b) Growth rate a of small harmonic perturbations vs 
the wavevector k for several values of the swapping ratio A (as 
labeled) in model M2. (c) Peak value of a vs swapping ratio 
A for model M2 (solid) and M3 (dashed). The growth-rate 
ratio Q is normalized by the characteristic time rj~^ — 7noK^, 
where D is the dimensionality of the model. 
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where Kg is the swapping range and Kc collision range. 
For a collision parameter A smaller than the swapping 
range Kg, particle swapping leads to the exchange of 
transverse particle positions, consistent with the results 
of our hydrodynamic calculations illustrated in Figs.[ljc) 
and[3l^a). Hence, the two terms of the integrand in Eq. 
([1]) cancel each other. For the collision parameter in the 
range Kg < |A| < Kc, particle interaction results in finite 
particle displacements. 

In MB, our 3D model [l3|, particles on coUisional tra- 
jectories behave as hard spheres without HI. In the swap- 
ping trajectory region |A| < Kg the particles exchange 
their transverse positions z (as in model M2). 

Small perturbation analysis - We show that the exis- 
tence of the swapping region | A| < Kg results in a layering 
instability, provided that the swapping ratio 



(a) r 



(3) 



is sufficiently large. The conditions for a uniform particle 
distribution uq to become unstable can be derived by 
analyzing small harmonic perturbations n(z;t) — hq + 
ni{t)e^'^^. For the translationally invariant population- 
balance models M2 and M3, such Fourier modes evolve 
exponentially, ni{t) = ni(0)e"*. 

FigurelSfb) shows that the growth rate a{k) has a peak 
around k w 2it/kc. Since the peak value becomes positive 
at a critical swapping ratio X — Xc, the uniform particle 
distribution is unstable to small perturbations for A > Ac, 
leading to formation of particle layers. Fig. [3Jb) depicts 
model M2, but we find that M3 yields a similar behavior. 

Figure [3l[c) shows the peak values of a plotted vs. the 
swapping ratio A (both for M2 and M3). For M2 the crit- 
ical swapping ratio is Ac = 0.24, and for M3 we obtain 
Ac — 0.213. Our BMC-HI simulations yield the instabil- 
ity in a similar parameter range. 

Based on the dynamics of pair collisions (with HI), the 
swapping ratio ([3]) can be controlled in two ways. First, 
the swapping range can be changed by varying the 
wall separation H. For example, for a particle pair in 
the middle of the channel, the dimensionless range Ks/d 
varies between 0.27 for H/d = 5 and 0.12 for H/d = 40 



11| . Second, the swapping ratio can also be controlled 
by changing particle roughness, because collision range 
Kc diminishes with the decreasing roughness amplitude 
|15| . as shown in Fig. EJa). We thus predict that by 
decreasing the magnitude of the particle roughness, we 
can induce formation of a long-range order in a confined 
suspension under shear. This prediction is confirmed by 
our BMC-HI simulations (cf.. Fig. [21 where A = 0.163 for 
e/d = 0.25 and A = 0.23 for e/d = 0.064, based on the 
mid-channel swapping range Kg/d = 0.189). 

Near-wall micro structure - Further predictions re- 
garding the wall-induced suspension ordering can be 
obtained when geometrical constraints that disallow 
particle-wall overlaps are introduced into the population- 
balance model (IlJ . Figures U and [S] show our results for 
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FIG. 4. Long-time near-wall suspension microstructure for 
subcritical values of swapping ratio A (model M3) (a) Density 
profile n/no for A = (dashed line) and A = 0.1 (solid), (b) 
inverse correlation length l^^ vs A. 



model M3 with such constraints. We consider a system 
with a large wall separation H/ Kc — 60 to emphasize 
two important regimes of suspension behavior: subcriti- 
cal (Fig.[4|) and supercritical (Fig. [5|). 

In both regimes we observe formation of a wall-induced 
layered order, but the evolution and extent of the layered 
microstructure is regime-specific. For subcritical values 
of the swapping ratio the microstructure appears 
only in the near-wall domain, whereas for supercritical 
values the layered structure propagates from the wall into 
the bulk of the suspension. 

Subcritical regime - Figure Slja) shows the subcriti- 
cal microstructure for A = (no swapping) and A = 0.1 
(significant swapping) . In both cases several particle lay- 
ers form near the wall due to the excluded- volume effect. 




FIG. 5. Near-wall suspension microstructure for supercritical 
values of swapping ratio A (model M3) . (a) Density profile at 
times ti/ro = 1250 (dashed line) and 12/ tq = 2500 (solid line) 
for swapping ratio A = 0.18 and channel width H/d = 60. 
Characteristic time tq — 1/ (no7Kj?). (b) Normalized propa- 
gation velocity Vf — VfTo/Kc- Near critical scaling behavior 
of w/ for two cases: (c) A > Ac, and (d) X> Xw 
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FIG. 6. Layer restructuring for supercritical swapping ra- 
tio A — 0.24 (model M3). (a) Merging of two central layers 
when the width is incommensurate H/kc — 15.4; interme- 
diate time (dashed line), long time (solid), (b) Number of 
layers vs normalized channel width: maximal number during 
evolution (solid); steady state (dashed); estimate based on 
the wavelength of the most unstable Fourier mode (dotted). 



The correlation length I of this local layered structure 
is insensitive to the value of A in the whole subcritical 
regime [cf., Fig. Hl^b)]. However, for A = the layers 
are not sharply defined, with many particles in the val- 
leys between the density peaks. For A = 0.1 the layers 
are much sharper. Thus the ST mechanism influences the 
form but not the extent of the subcritical microstructure. 

Supercritical regime - In the supercritical regime (cf. 
Fig. [5]) the wall-induced microstructure propagates from 
the wall with velocity w/, forming a growing number of 
particle layers, as shown in Fig.[SJa). In a fully developed 
microstructure the layers are well separated, and particle 
jumps between layers are rare [cf., Fig. [IJd)]. 

The propagation velocity i"/ of the layered structure 
in this regime [cf. Fig. [SJb)] is an increasing function 
of A, with two distinct dynamical domains. For A > Ac 
(suspension unstable to small perturbations), Vf rapidly 
increases with the swapping ratio, showing a power law 
critical behavior Vj ^ (A — X,^^!"^ [cf. Fig. [Sfc)]. In 
the domain A^, < A < Ac the suspension is stable to 
small perturbations, but unstable to the large perturba- 
tion caused by the wall. In this regime, the data can 
be well fitted to the power law vj ^ (A — A^,)'^ with 
A^ w 0.11 and /3 « 4 [cf. Fig. M^d)]. Due to the large 
value of the exponent /3, the propagation velocity w/ is 
practically zero for A < 0.16. 

The layer separation in the propagating microstruc- 
tural fronts may be incommensurate with the wall sepa- 
ration. In such cases the two mid-channel layers merge, 
as depicted in Fig. IHl^a). After such suspension restruc- 
turing, the width of the density peaks significantly de- 
creases. Note that the initial number of layers agrees 
fairly well with our prediction based on the wavelength 
of the most unstable Fourier mode [cf. Fig. Hl^b)]. 

Conclusions - The ST mechanism significantly alters 
the suspension microstructure in a confined Couette flow 
by induciiig layering. Previously (solving a long-standing 
paradox [16||). we have demonstrated that swapping tra- 



jectories cause the anomalous enhancement of hydrody- 
namic diffusion, and that they stabilize particle chains in 
microfluidic channels 1]|. The layering behavior iden- 
tified and analyzed here adds to the growing evidence 
that a seemingly subtle ST effect has far-reaching conse- 
quences in confined particulate flows. 

Our direct numerical simulations indicate that swap- 
ping trajectories significantly influence a variety of con- 
fined systems, including suspensions at higher concentra- 
tions [IJI and flow-driven particle monolayers described 
in Q- Since such monolayers (in the disordered state) 
have a similar structure to the layers observed experi- 
mentally in [l^, we conclude that the ST effect has sig- 
nificant implications for suspensio n dy namics. In suspen- 
sions with nonzero inertial forces [17| and in viscoelastic 
fluids, swapping (reversing) trajectories occur even with- 
out confinement, so the ST mechanism is likely to affect 
suspension microstructure also in such flows. 

Finally, we note that the suspension behavior described 
in our study significantly differs from the dynamics of a 
dilute gas. In a gas, binary collisions always produce 
a homogeneous equilibrium state, according to Boltz- 
mann's H-theorem. By analogy, it is assumed that binary 
collisions in suspension flows always lead to a diffusive be- 
havior that results in relaxation of density fluctuations. 
Here we report that uncorrelated binary collisions pro- 
duce inhomogeneous layered microstructure. 
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